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Abstract 

For a dissipative system with Ohmic friction, we obtain a simple and exact solution for 
the wave function of the system plus the bath. It is described by the direct product in two 
independent Hilbert space. One of them is described by an effective Hamiltonian, the other 
represents the effect of the bath, i.e., the Brownian motion, thus clarifying the structure of the 
wave function of the system whose energy is dissipated by its interaction with the bath. No 
path integral technology is needed in this treatment. The derivation of the Weisskopf-Wigner 
line width theory follows easily. 

The simplest example of a dissipative system, an harmonic oscillator coupled to the environ- 
ment, which is a bath of harmonic oscillators, has been the subject of extensive studies[l- 15] . 
We shall show in the present paper that in a special case, the Ohmic case (to be defined later), 
the dissipative system can be exactly treated both classically and quantum mechanically, thereby 
clarifying the sense in which the wave function is describable by an effective Hamiltonian. In this 
treatment path integral technology is not needed, and our presentation is self-contained. 

We consider the problem discussed by Caldeira and Leggett (CL)[l], an harmonic oscillator 
system (the dissipative system) with coordinate q, mass M, and frequency ( uf^ + Acj 2 ) 1 / 2 , inter- 
acting with a bath of N harmonic oscillators of coordinates Xj, mass m^, and frequency u>j., where 
Au> 2 is a shift induced by the coupling already discussed by CL. The Hamiltonian of the system 
and the bath is: 

H = 2M + \ M ^ 2 ° + Auj2 ) q2 + + 5 (!) 

The dynamic equation for operators in the Heisenberg representation leads to the following set 
of equations of motions: 

Mq = — MAu 2 q — ^ CjXj , ( 2 ) 

3 

rrijX = —mjbJ 2 Xj — c 3 q ( j = 1,2 ( 3 ) 
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Now, applying the Laplace transform[2] (the bars are our notations for the Laplace transform, s 
is the Laplace transform of time t), equations (2), (3) can be used to eliminate bath variables Xj 
to obtain the equation for q 


M(s’t - sqo - ») = —Mtiilq - MAu’ij - £ - E 


9, 


(4) 


where qo , qo, x j0 , x j0 are the initial values of the respective operators in the Heisenberg represen- 
tation. Assuming the number of bath oscillators is large enough so that we can replace the sum 
over j by an integration over uj, the coefficient of the last term can then be separated into two 
terms: 


f°° p{wj)cjduj _ 2 r 
Jo rrij u 2 Jo 


00 pMtfdujj 


- m : u](s 2 + w 2 ) ’ ^ 

where p(u>j) is the bath oscillator density. Following an argument similar to the one pointed 
out by CL[1], the requirement that the system becomes a damped oscillator with frequency u Q 
and damping rate q in the classical limit, known as the “Ohmic friction” condition, leads to the 
following constraint: 

2i ]M rrijoj] 


PM = 


7T 


( 6 ) 


By observing eq.(4) and eq.(5), it can be shown that the second term of eq.(5) leads to damping 
with the damping constant rj , while the first term of eq.(5) represents a frequency shift. If the 
frequency renormalization constant A uj 2 is chosen to satisfy: 


MAu J 


f°° pMcjdujj 


-L 


rrijLO 2 


(7) 


the frequency is shifted to u>o- Then equation (4) is simplified, and its inverse Laplace transform 
gives the quantum Langevin equation, valid at time t >= 0+: 


q + qq + u 0 q 

with the Brownian motion driving force: 

m 


m , 


• sin uJjt 

- 2^ Cj(Xj 0 cos Ujt + Xjo ) . 


( 8 ) 


(9) 


During the derivation, in order to carry out the integral in eq.(5), we used the requirement of the 
inverse Laplace transform that s must pass all the singular points from right of the complex plane, 
and hence real(s) >0. 

Equation (8) and (9) are the equations of a driven damped harmonic oscillator, the solution 
of which is well known as a linear combination of the initial values at 90, < 7 o, Xjo, and x j0 : 

q(t) = ai(t)q 0 + a 2 (t)q 0 + Y, + b j2 (t)x j0 ) , (10) 

j 

Xi(t) = an(t)q 0 + a i2 (t)qo + ^(/?,ji(<)^jo + Pij 2 (t)ij 0 ) , with (11) 

j 
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( 12 ) 


= 


ve 


-vt 


fie 

fi 


— ut 


■» a 2{i) = 


e _Mt - e -l/t 

V — fl 


ViV- 
fi = — + u> i , and v — on , 

2 2 


(13) 


and here u> = (o>q — is the frequency shifted by damping. (All formulae are correct whether 

u is real or imaginary. To avoid a minor detail of the initial value problem, we have redefined the 
initial time as t=0+.) The explicit expression for bji, bj 2 , a,i, a; 2 , Aji> and /?, y 2 are well known 
in freshmen physics. 

We emphasize that the use of the Laplace transform instead of the Fourier transform allows 
us to express q(t) and Xj(t) explicitly in terms of the initial values, as in eq.(10) and eq.(ll). 

The equations (10) and (11) serve as the starting point of subsequent discussions. We will 
proceed to find the Green’s function of the full system, and hence the solution of the wave function 
in Schoedinger representation. The result tells us in what sense the damped oscillator is described 
by an effective Hamiltonian without the bath variables and gives its specific form, it also shows 
that under this condition, the wave function can be factorized, and the main factor relevant to 
the damped oscillator is a solution of the Schoedinger equation with an effective Hamiltonian. 

Equations (10) and (11) are correct both in classical mechanics and in quantum mechanics in 
the Heisenberg representation. We notice that q(t) and ay(t) are both linear superpositions of 
? 0 i 9o, Xjo, ijo with c-number coefficients. The commutation rules between q(t), q(t), Xj(t), xj(t) 
are = m 1 ^i(^)] = ”■> and operators q(t) and q(t) commute with Xj(t) and 

ij(t). One can prove these commutation rules by two ways: (a). By direct computation, using 
the fact that at t=0, they are correct, (b). By the general principle that q(t), q(t), Xj(t), ij(t) 
are related by a unitary transformation to q 0 , q 0 , xj 0 , andijo. 

Equation (10) and (11) show that the operators q(t) and Xj{t ) can each be written as a sum 
of two terms: 

= + (14) 

3 


x i(l) — (i(t) T ]E ^0(0 > (15) 

J 

where Q(t) and £(£) are linear in q 0 and q 0 and independent of Xjo and ij 0 , and and X,j(t ) 
are linear in x j0 and x j0 , and independent of q 0 and q 0 . Thus Q(t) and (i(t) are operators in 
one Hilbert space Sq, while £j(t) and Xij(t) are in an independent Hilbert space Sx, and the full 
Hilbert space is a direct product Sq <g> Sx- 

We shall first analyze the structure of Sq. We write that: 


Q(t) — aiQo + a^Qo — Q-iQo — 0-2 


ih d 

MdQo ‘ 


(16) 


To explicitly show that we are discussing the Sq space, we define Q 0 = q 0 , Q 0 = q 0 . The eigen- 
function of Q(t) with an eigenvalue denoted by Q x , in the Q 0 representation, is easily calculated 
to be: 


u Qi(Qo,t ) 


2nh sintuf 


exp 


' iM 
2ha2 


(aiQl-2Q 0 Qi+4>(Qi,tj) 


(17) 
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with <f> as an arbitrary phase, i.e., a real number. This eigenfunction is related to the Green’s func- 

tion G{Q u Qo\t,0) =< Qi\U q (t)\Q 0 >, by u Ql {Q 0 ,t) =< Qol^WIQi > =< Go|tf£(0IGi >= 
G*(Qi, Q 0 ] t, 0), where we denote the evolution operator by U(t), which is unitary when we choose 
the eigenvectors of Q(t) to be orthonormal. Thus we have: 

G{Qi,Qo;t,0) = exp[— (a 1 gg-2Q 0 C?i + Wi,0)] • (18) 

Next, we shall determine the arbitrary phase <f>{Qi,t), which is the phase of the eigenvectors 
of Q(t). Using eq.(16), we find the commutation rule for Q and Q : 

e~ nt 

[Q,g] = ( fl id2 — a ia2 ){Q 0 ,Q 0 } = e ^[goi^o] = • (1^) 

Thus we define the canonical momentum P(t) as: 

xft d 

P{t) = Me vt Q(t) = M e^dxQo + d 2 g 0 ) = diQ 0 - > (20) 

and get the commutation rule: [Q(t),P(t)] = th. The eigenfunction of P(t) can be calculated in 
two ways: (a). We can calculate the eigenvector of P(t) in the Qo representation using eq. ( 20 ) 
and then use the Green’s function eq.(18) to transform it into the Q(t) representation; (b). The 
commutation rule [Q(i),P(t)] = ih requires that P(t)=-*7iJ^, so the eigenfunction of P(t) with 
eigenvalue Pi is exp[i^fQ]. By comparing these two solutions, the arbitrary phase 4>{Q i,t) in the 
Green’s function is determined to within a phase <j>(t), which is independent of Q\. is an 
arbitrary real function of time, except that 0(0) = 0 so that it satisfies the condition that at t=0, 
the Green’s function becomes 6(Qi - Qo)- Thus we obtain the Green’s function in the Sq space: 


G{Qii Qo ; t, 0) — 


Mues’’* 
2irih sin ut 


exp 


2 ha 2 


(a 1 Q 2 + h 2 e r}t Q 2 1 -2Q 0 Qi) - 




( 21 ) 


It is then straight forward to derive the Hamiltonian Hq using the following relation: 

Hq = , ( 22 ) 

and remembering that the matrix elements of Uq and U q 1 are the Green’s function and its 
conjugate. The result is: 

H Q = e~ , ’ t ^+ l - M u; 2 0 e'' t Q 2 + m . (23) 

Since 0 is arbitrary except that 0(0) = 0, we can take <f>(t) = 0. Therefore we have derived the 
well known effective Hamiltonian for the dissipative system. We emphasize that the expression for 
Hq is here derived, while in usual literature it is introduced by more or less heuristic arguments. 

Next, we shall analyze the effect of the bath. Similar to eq. (17) we obtain the eigenfunctions 
0t }i ( x jo,t) for (j. Using Dirac’s notation we have: Q\u q ^ >= >, >= >• 
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Thus |uq > 1 0^ > is an eigenvector of q(t), with eigenvalue of Q + ^ In other words, 

•»' i 

the eigenvector of q(t) with eigenvalue q is \q, {£,} >= |u _ ^ > II 0 1^ >• 

i 

If initially the wave function is I'I'o >= IV’o > H (^) |xjo >, to calculate the wave function at 

3 

time t, we should expand in terms of the eigenvectors |<f, {£,} >, i.e., we should calculate 

=< uqIV’o >= J UQ(q 0 )ipo(qo)dqo , (24) 

Xi(£i>^) =< '- I Xjo > J (3'jo)Xjo(&jo)dxjo . (25) 

Then the wave function in the Schoedinger representation at time t is 

nq,{tj},t) =< q, {&} |*o >= V>(9-£ £,-,*) UxAM • (26) 

j j 


Notice that of eq.(24) is the wave function in the Schoedinger representation with the 

effective Hamiltonian eq.(23). Hence we have connected the effective Hamiltonian approach to the 
dissipative system problem with the other approaches that take both the system and the bath into 
account. We also notice that even though our $( q , {£j}, <) is in a different representation from that 

°f the usual probability interpretation is still valid: J J ... J |ip(?, {£/},<)| 2 is 

the probability density of finding the particle at q. Since this solution is very simple, it provides a 


simple way to analyze other complicated problems, e.g., study the influence of Brownian motion 
on interference, which we shall not elaborate because of space. 

Under certain conditions, eg., at low temperature and when the system q is in highly excited 


states so the range of q is large enough that for all the values of which do not have vanishingly 
small probability, q > 1 we can approximately write: \f(g, {£,*},<) = 0 0* 


That is, the wave function is factorized, the dissipative system q can be described by the wave 
function on ly> and the Brownian motion can be ignored. Therefore it is interesting to 

examine the width of the argument of the wave function 0, due to the Brownian motion, i.e., 
the mean value of (^£?) 2 at time t. It can be calculated using its expression in the Heisenberg 
j 

representation as introduced by (10) and (18). At temperature T, this width is : 


<(EfeW) 2 >=E^(feWI 2 + ^M<)| J )«><(^) . (27) 

This width is zero initially, then approaches its final equilibrium value in a time interval of the 
order of !/»/. At low temperature limit, the equilibrium width is simplified to: 


a 2 {t = oo) = 


h 

2xmu> 


— + arctan 

4 



(28) 


If the damping rate rj is much smaller than the frequency of the oscillator u> 0 , this width happens 
to become the same as the width of the ground state of the system h/(2mu). 
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Finally, it is interesting to see the distribution of the dissipated energy of the harmonic oscillator 

in the bath, and check if it agrees with the Weisskopf-Wigner line width theory[16] For simplicity, 

we assume zero temperature. To calculate the energy dissipated by the system into the j’th bath 

oscillator, we use eq.(ll) and its derivative to obtain the expression of Xj(t) and pj = nijXj^), 

p' 2 

which are then substituted into the expression of the energy of the j’th oscillator: hj = + 

IrrijUJjXj. We then calculate the expectation value of hj, assuming initially the system is in the 
n’th excited state, and the bath oscillators in the ground state. We calculate the contribution to 
the expectation value of hj from the system, by keeping only terms which depend on q 0 , and q 0 . 
The result is then the energy disspation by the system into the j’th oscillator. When multiplied 
by the density of states p(<x>j) eq.(6), it gives the dissipated energy spectrum. It is a function of 
Uj, with a narrow peak near the resonance u>j = u, if the damping is small, i.e., if r) u. It is 
oscillatory with a frequency of 2u>j. Its time average over a period is : 


_ . 1.. n 1 

Ej = (n + -)hu> —— 

2 <r£ + ( w - i »;) 2 




(29) 


where A(u>j) varies slowly near the resonance, i.e., over the width t// 2 of the peak it changes very 
little. Therefore, we can replace Uj by u> near the peak, and the result is simplified to: A(ui) « 1. 
Thus eq.(29) shows that the dissipated energy has a Lorenzian distribution near the resonance, in 
agreement with the Weisskopf-Wigner line width theory. 
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